A random matrix theory of decoherence 

T Goring C Pineda^-^S H Kohler^ and T H Seligman^^ 

^Z. • ^ Departamento de Fisica, Universidad de Guadalajara, Blvd. Marcelino Garcia 

f~^ I Barragan y Calzada Olimpica, Guadalajara C.P. 44840, Jalisco, Mexico. 

Cs| ■ ^ Institut fiir Physik und Astronomie, University of Potsdam, 14476 Potsdam, 

Germany 

•^ Instituto de Ciencias Fisicas, Universidad Nacional Autonoma de Mexico, Mexico 
^— V ■ ** Centro Internacional de Ciencias, Cuernavaca, Mexico 

ff^ I ^ Fachbereich Physik, Universitat Duisburg-Essen, D-47057 Duisburg, Germany 

E-mail: thomas.gorin@red.cucei.udg.mx, carlospgmat03@gmail.com 

Abstract. Random matrix theory is used to represent generic loss of coherence of 
Jlj ' a fixed central system coupled to a quantum- chaotic environment, represented by a 

2 . random matrix ensemble, via random interactions. We study the average density 

^H' matrix arising from the ensemble induced, in contrast to previous studies where the 

average values of purity, concurrence, and entropy were considered; we further discuss 

when one or the other approach is relevant. The two approaches agree in the limit of 

^ ' large environments. Analytic results for the average density matrix and its purity are 

CO , presented in linear response approximation. The two-qubit system is analysed, mainly 

—J. I numerically, in more detail 
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1. Introduction 

Decoherence has become a subject of increasing interest, since it is not only related 
to the problems of measurement and emergence of classical behaviour [H [21 El H] 
but has become one of the main stumbling blocks for quantum information solutions. 
Under most usual circumstances decay of coherence is exponential and this has been 
experimentally confirmed by experiment [5], [6] . Exponential decay of coherence typically 
holds if the Heisenberg time of the environment is to long to provide a relevant time 
scale, i.e. if decoherence is complete long before the Heisenberg time. As the global 
environment is always present and has, for all practical purposes, infinite Heisenberg 
time, our discussion seems spurious but this is not the case. Quantum information 
tasks require sufficiently low temperatures, high vacuum, good screening and careful 
selection of the states used etc, such that the global environment can be neglected 
to first approximation. We are left with a near environment consisting of only a few 
degrees of freedom which are involved in the control of the operations of the quantum 
gates. This near environment may have very specific properties, but is often governed 
by dynamics that might go for "quantum chaotic". In this paper we provide a generic 
description of decoherence which also covers such a situation. 

This requires modelling of the environment, while avoiding a complete microscopic 
description. More specifically, we describe a generic environment with "quantum 
chaotic" dynamics by random matrix theory (RMT). The first step in this direction 
is the pioneering paper of Lutz and Weidenmueller [TJ E] in which they propose a RMT 
description that reproduces the results of a Caldeira-Leggett model. Thus they show 
that RMT also works in the standard situation. We shall see that this is closely related to 
choosing a very dense spectrum in the environment, which sends the Heisenberg time to 
infinity, or, more generally, where decoherence is strong, long before the Heisenberg time 
of the environment is reached. While this is very reasonable for the far environment, the 
near environment can have level densities comparable to those of a central system, and 
this is where RMT models can produce new results. A model to study such effects 
was first proposed in [9l [10] and further developed in [HI [12]. These models are 
closely related to one developed for fidelity decay [I3], which was successfully tested 
in experiments [HI, [151 [IB]- For decoherence, the models were compared to numerics 
with a spin chain environment pTT] . 

In these models purity was calculated analytically in linear response approximation 
as averaged over the random matrices. For the purpose of quantum information, this 
was quite sufficient, as only high purities are needed. Yet, the beauty is, that these RMT 
models are sometimes exactly solvable by super-symmetric techniques. An attractive 
example are the solutions for the decay of the fidelity amplitude [HI [TH [201 [21] . 
Yet, these techniques nowadays are largely limited to the calculation of two-point 
functions. Purity is obviously of higher order, and thus we should look for an amplitude. 
Coherences, i.e. off-diagonal elements of the density matrix are good candidates, and 
their decay is often used to describe decoherence [5|. This suggests to consider the 
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average density matrix, something that has not been done so far. 

Actually, in the entire formulation of the RMT models a basic question has not 
been addressed. If we consider the evolution of some initial state by an ensemble 
of Hamiltonians we also get an ensemble of density matrices for the central system. 
Normally, we would consider the average density matrix as the basic quantity from 
which observables such as purity or concurrence should be calculated. However, since 
these quantities are non-linear in the density matrix, their average (over the ensemble 
of density matrices) is different from their value computed from the average density 
matrix. At the end, it is the particular application which determines which kind of 
averaging must be chosen. We shall discuss this in section [5], but we will find in the 
previous sections, that for large environments the quantities we analyse, namely purity, 
entropy and concurrence converge to the same result for both averages. 

We first formulate the random matrix model we use and derive a few basic relations, 
which will be needed in subsequent sections. Some obvious conditions the model must 
fulfil are presented before proceeding to calculate the average density matrix in linear 
response approximation. For the particular case of one qubit as a central system we 
outline how to obtain the exact answer with super-symmetric theory. Similar results for 
the average density matrix are presented in what we earlier called spectator configuration 
[22] . This configuration has proven very useful to describe more complex central systems 
with high purity [23]. We next concentrate on two-qubit central systems in the spectator 
configuration, where we display some properties of the density matrix ensemble including 
statistics of the eigenvalues and some further results on a curious relation to Werner 
states shown in previous papers [12l [T7]. In section [5] we discuss in more details the 
meaning of an ensemble of density matrices. The discussion will include both situations 
where one should use the average density matrix to calculate non-linear functionals and 
others where this is not appropriate. This will lead us to conclusions and an outlook on 
the relevance of more realistic RMT environments. 

2. Random matrix models 

We shall start with a central system that does not display any particular structure, and 
thus the entire system will be bi-partite consisting simply of this central system and an 
environment. Decoherence will be viewed as the entanglement with the environment 
after unitary evolution of the entire system. As we wish to describe a near environment 
with quantum chaotic features, we neglect the far environment entirely, as indeed 
we shall throughout this paper, assuming that decoherence resulting from the far 
environment occurs on a much longer time scale. 

The Hilbert space is divided into two parts, that of the central system Tic and that of 
the environment Tie- The dimension of those Hilbert spaces is denoted hj m = dim(Hc) 
and A^ = dim(7ie), respectively. Eventually, the case N ^ oo will be particularly 
important. For a factorized initial state of the form 

^0 = ^c ® ^e (1) 
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we solve the Heisenberg equation 

idtg{t) = Hxg{t)-g{t)Hx, (2) 

where the Hamiltonian is of the form 

Hx = H, + H, + XV , (3) 

such that Hq = H^ + H^. Note that typically, we shall not require the initial state to 
be pure. It may well be a product of two mixed states ^c and Qe- In the following 
step, we express the forward evolution of the whole system in terms of echo dynamics. 
Echo dynamics means, that we consider a forward time evolution which includes the 
coupling between central system and environment and a backward evolution, which is 
driven by the internal dynamics of environment and central system alone. The latter 
will not modify the entanglement between the two. This is important, because then we 
may apply the linear response approximation to the echo-dynamics which results in a 
much larger range of validity. A systematic application of the interaction picture would 
do the same simplification. Hence, we write: 

0{t) = e-^^* go e'^^* , (4) 

where we have set h = 1. In order to obtain the echo dynamics, we evolve the system 
backwards in time with the Hamiltonian Ho: 

QMit) = M,(t) ^0 M,(t)t M,(t) = e^^°* e-^^* . (5) 

This may be interpreted as a forward evolution with the echo-operator M\{t) ^T]. The 
evolution of the central system is then given by 

Qc{t) = tre [ g{t) ] = tre [ {u^ ® Me) QAlit) (mJ » «!) ] , (6) 

where Uc = exp{—iHct) and Ue = exp{—iHct) . This expression for ^c(^) rnay be further 
simplified: 

gc(t) = Uc tre [ (Ic (S) Me) QM{t) (Ic ® ul) ]ul = Uc Qcit) mJ , (7) 

where ^c(^) = tYc[gM{t)] is the density matrix of the central system in the interaction 
picture. Note that for operators of the form Ic ® Me the cyclic permutations are allowed 
even though we are dealing with a partial trace. Equation (jTj) holds only as long as Ho 
describes uncoupled dynamics. 

Linear response approximation In the linear response approximation, the echo-operator 
Mx{t) reads 

Mx{t) ^ 1 - iA I{t) - AV(t) , (8) 

where 

I{t) = f dr V{t) , J{t) = f dr f dr' V{t) ^(r') , (9) 

Jo Jo Jo 

and V(t) = exp{iHot) V exp{—iHot). Starting from equation ([8]), we obtain the density 

matrix of the central system averaged over the coupling matrix V (this average will be 
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denoted by angular brackets). Subsequently, we may also average over the environmental 
Hamiltonian H^. However, for the linear response result as such, this is not necessary. 
In the interaction picture, we therefore obtain: 

(g^t)) ^ tre ( ( U - iA I{t) - \^J{t) ) f?o ( 1 + iA I{t) - \^J{t)^ ) ) 

^Q,- \\ {Aj) - {A,) ) , (10) 

where we have used the notation 

Aj = tl,[J{t)QQ + QQj{t)^], Ai = il,[I{t)QQl{t)]. (11) 

Note that I{t) is self-adjoint, while J{t) is not. According to Equation (jTj), this yields 
for the time evolution of the central system: 

(f?c(t))=Mc(t) (^c(t))Mc(t)^ (12) 

The averaging procedure itself as well as the different possible statistical assumptions, 
will be discussed below. 

Purity We can now obtain the purity of the average density matrix of the central 
system. In view of ([7]), we find for the purity of the averaged state ^c(^): 

P(t) = tr (f?e(t))' = tr[Me (Ut)) uluc (^c(t)) mJ ] = tr (^e(t))' . (13) 

The purity may be more explicitly written in terms of the average state of the full system 
(in the interaction picture) as 

Pit) = tre ( tre (^A/(t)) tre (f?A/(t)) ) . (14) 

Note that in this expression, we cannot exchange the order of the partial traces, since 
{0Ai(t)) is no longer a pure state in the product Hilbert space H of central system and 
environment. This is the reason why we cannot simply adopt the calculations in [22] to 
the present case. In linear response approximation we find with the help of (TTO!) 

P{t) ^ tr( qD-X^ tr( g, (Aj) + (Aj) g, - g, (Ai) - (Ai) g, ) 

^ tr( gl ) - 2X' [tT{{Aj) g,) - tr((A,) g,) ] . (15) 

Since Aj, Aj and ^c are all Hermitian operators, these expectation values are real. 

Note the difference to the approach chosen in Refs. P, [171 [22l [23], where we 
computed the average purity rather then the purity of the average density matrix. This 
difference will be discussed as we proceed and taken up in detail in section [5l 

Statistical (random matrix) assumptions We divide the statistical assumptions 
employed into two parts. A basic part which is invoked when we compute the density 
matrix of the central system, and an optional part, which makes stronger assumptions 
on the dynamics of the environment in order to obtain more explicit results. 

• We assume that V is taken from the Gaussian unitary ensemble. This allows to 
work in the eigenbasis of Hq, where Hq is diagonal. The transformation into this 
basis is necessarily a unitary transformation, such that the ensemble for V remains 
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unchanged. In addition, due to the separabihty of Hq = H^ + He, the eigenbasis 
may be chosen as the product basis \ij) of eigenstates of H^ and H^, respectively. 

• The choice of H^ is quite arbitrary as long as it remains fixed. In contrast, we 
eventually assume H^ to be randomly chosen from an appropriate ensemble, say 
one of the classical ensembles. In that case, the two-point form factor 62 (^) [21] is 
the only characteristic quantity of the spectrum. The two-point form factor may 
well be left unspecified. However, as specific examples we may consider Poisson 
spectra (62 (^) = 0), as well as GOE- or GUE-spectra [21] . 

2.1. Average density matrix of the central state 

Here, we calculate the behaviour of the central system, when averaging over the random- 
matrix coupling to the environment. The state of the central system is then described 
by the density matrix {gdt)). In view of Eqs. ( TTOl) and ( TT2l) . we start by calculating the 
terms (Aj) and (Aj). 

{Aj) = f dr r dr' tre [ {V{t) V{t')) q, + ^0 {V{r') V{t)) ] , (16) 



where 

{V{t)V{t')) = \zj) e'(^---^"--') (F.,,HVfc,,„„)e-'^-(---') {mn\ 

= lu) ^e-^(^^'-^'^)("-"') (zjI = C,{t - r') ® C^ir - r') . (17) 

kl 

We use the implicit summation convention for indices appearing more than once. The 
diagonal matrices Cc and Ce are defined in analogy to Ref. [22] with the only difference 
that we use the subscript c instead of 1 to refer to the central system: 

air) = J2 K) e-'(^-^>)- {^\ Ce{r) = J^ \j) e-'^^'"^^)- {j\ . (18) 

ik jl 

The separable result in fITTj) is due to the choice of the GUE as random matrix ensemble 
for the coupling. In the GOE case, we would obtain an additional term, which breaks 
the separability. 

{Aj)= [ dr rdr'[C7e(r-r')Ce(r-r')& + herm.conjg.], (19) 

Jo Jo 

where Cx{t) = tTx[Cx{T) g^] for a; = c,e, again in full analogy to Ref. [22j. For the 

second term {Aj), we obtain 

{Aj) = [[ drdr' tre [ (^(r) go V{t')) ] , (20) 



where 

tre [ (Vir) go V^(r')) ] = \t) e'(^»^-^«)- {V,,,m V^,,,) gl^ gl e^(^™-^-)-' {p\ 
= K) e^(^-^-^-)- gl^ gl (5,/,^5,„ ^^iE^r.-E,,y ^p\ 
= \i) e^(^'^-^«)- gl^ gl ei(^'='-^'^)"' (i| 
= a(r'-r)|.)e-(^-^«)(— ')^^,(z|. (21) 
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This expression is a weighted version of Cc(t), with weights given by the diagonal 
elements of Qc- It may be understood as a mapping of the density matrix g^ onto a 
diagonal matrix. It will be denoted by 

C,(t) : ^e ^ CMit) = E 1^) e-'(^^-^^)* gl, {^\ • (22) 

ik 

We will encounter a generalization of Cc{t) again in the following section. Note that 

trCc[^c](t) = Cc(-t). (23) 

We may write: 

{Ai) = jj drdr' Ce(r' - r) C,[Q,]{r - r') . (24) 

Note that the double integral (Aj) is Hermitian though the integrand is not. We thus 
have obtained the relevant quantities for computing the average density matrix in terms 
of correlation functions. In order to obtain an explicit expression, we can subsequently 
average over the environmental Hamiltonian also. 

2.2. Master equation 

In order to arrive at the desired master equation, we will make two approximations. 
First, we assume that the time scale for decoherence is much shorter than the Heisenberg 
time th of the environment. We then approximate the spectral correlation function in 
the environment with a simple delta function: 

Ce(t) -. Sit/TH) . (25) 

In the case of the random matrix environment Ce{t) — > Ce{t), the latter also contains 
a smooth part, which becomes important at times of the order of th and beyond; 
for the precise form of C'e(t) see j22]. Second, we assume that the state g(t) of the 
full system (central system and environment) remains a product state throughout the 
whole evolution. One common argument is that Qc does not change in the course of time 
because the environment is much bigger than the central system [25]. Together, these 
approximations are consistent with a Markov approximation. We then use the linear 
response result to advance the density matrix {gdt)) for a time step, which should be 
small compared to the decoherence time but large compared to the mixing time in the 
environment. We thereby obtain a closed differential equation for {gdt)), the desired 
master equation. 

As before, we restrict the discussion to the GUE coupling. Note though that a 
GOE coupling can be treated along the same lines and would yield precisely the same 
result. In this respect, on the level of a master equation approach, the evolution of 
the central system does not distinguish between systems, where the coupling breaks the 
time-reversal symmetry and others that do not. This is a general feature of the Fermi 
golden rule approximation. With ( l25l) and due to Cc(0) = m Uc, we find: 



{Aj) = dr / dr' 5 i j [ Cc(r - r) ^c + herm. conjg. ' 
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= ^ dx[ Cc(0) Qc + herm. conjg. ] = m t Th Qc , (26) 

(Aj) = f drf dr'S (^^) Cc[&](r - r') = r^ f ^" dx Ue = r^ t Uc , (27) 

where we have used the fact that Cc[£ic](0) = Ic- With these two results and ( TTOl) . we 
obtain for the derivative of (^c(^)) at t = 0: 

^ (^c(t)) = -A' (mr^fPc-THlc) . (28) 

Since we are working in the vicinity of t = 0, we may replace Qc with (^c(O)) on the 
right hand side and thereby close the above relation: 

^ (^c(t)) = -m TH \' ( iQcit)) - m-' 1, ) . (29) 

We therefore obtain as a closed evolution equation for the density matrix of the central 

system: 

^ (f?,(t)) = -i [H,, {g,{t))] -ttithX^ (&(t)) - m-i U, ) , (30) 

which constitutes the desired master equation [25| . 

2. 3. Purity of the average density matrix 

We calculate the purity on the basis of (TT5|) . That implies the calculation of two terms, 
tr((Aj)pc) and tr((A/)^c)- For the first term we obtain from fll9j) : 

tT{{Aj)g,) = I dr I dr' tr[Ce(r - r') C,{t - r') qI + herm. conjg. ] 
Jo Jo 

drdr' Ce(r - r') tr[ Cc(r - r') gl] . (31) 

'0 

If Qc is a pure state, qI = Qc and 

ti{{Aj)Q,) = jj drdr' Ce(r - r') C,{t - r') . (32) 

Similarly, we obtain: 

ti{{Ai)g,) = jj drdr' C,{r - r') tr[Cc[f?c](r - r') f?,] 
t 



drdr' Ce(r' - r) S,{r - r') , (33) 

J Jo 

where 

^c(t) = $^e-(^-^-)*^^,^?,. (34) 

If Qc is a pure state, the quantity S'c(r — r') agrees with that defined in Ref. [22]. Overall, 
we obtain 

P{t) ^l-2\^ [ dr [ dr' [Ce(r - r') Cc(r - r') - Ce(r' - r) 5c(r - r') ] . (35) 

Jo Jo 
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For a large environment (A^ -^ oo) and for a pure state Qc, this result agrees with the 
average purity calculated in Ref. [22]. For finite N, their difference is calculated in 
Appendix A. Note that neither result requires any averaging over the dynamics in the 
environment. The averaging over the coupling is sufficient. 

2.4- Towards exact expressions for the averaged Density matrix 

For small dimensions of the central system exact expressions for the ensemble averaged 
entries of the reduced density matrix can in principle be obtained by methods similar to 
the ones used for the calculation of fidelity decay [201 EB [26] . For simplicity lets assume 
the central system to be a two-level system (i.e. one qubit), initially pure. We write 

3 

Pc(t) = ^Pcn{t)o-n , (36) 

n=0 

where cxj, i = 1,2,3 are Pauli-matrices and ctq is the unit matrix. Then the ensemble 
average of pcnit) takes the form 



p,nit) = tr [e-i^At(p^ (g) p^) ei^At(a„ ® 1^)J . (37) 

We use an over-line to denote the ensemble average, because here, the average is not 
only over the coupling but also over the environmental Hamiltonian and even over initial 
states. Due to the latter, the above expression further simplifies to 

3 ^ 

Pn{t) = Y,Pcra{^)fnm{t), /n„^(t) = ^ tr [c-i^A* g-'-^A^"*] . (38) 

m=0 

The quantities fnmit) are formally average fidelity amplitudes or -equivalent ly- the 
ensemble averages of the trace of echo operators. The forward time evolution is given 
by the 2A^ x 2A^ matrix amH\an and the backward time evolution is propagated by Hx. 
Fidelity amplitudes have been calculated exactly for several different random matrix 
ensembles using supersymmetry. The calculation of fnm{t) is therefore amenable to an 
exact treatment using the techniques of [211 [23 [2B] • We shall present the result in a 
future paper. 

3. Decoherence within the spectator model 

Dynamical model In the spectator model, the central system is divided into two parts 
Tic = 'Hi® 0.2) where only one part is coupled to the environment (in our case this will 
be Til). The Hilbert space of the environment is as before Tie. The dimension of the 
different spaces are: rrii = dim(7ii), m2 = dim(7i2), and A^ = dim(7ie)- The dimension 
of the total Hilbert space is A^m, where m = mi m2- The Hamiltonian reads 

H = Hi+H2 + He + Vi,e ® I2 . (39) 

For an initial state Qc ® ^e this means only Hi and He are coupled dynamically, whereas 
H2 may affect the dynamics of the whole system only due to an initial entanglement, 
i.e. anon-separable Qc- 
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Linear response approximation The spectator model may be treated along the same 
lines as the non-structured model considered in section [2l Thus we may use Eqs. fITOl) 
and (IT2|) without change, such that the principal task is to compute (Aj) and (Aj) on 
the basis of (fT6|) and (120|) . respectively. Hence, it is only the coupling matrix in the 
interaction picture together with its second moment, what we need to calculate anew. 
To that end we introduce Ui(t) = exp(— iifit), U2(t) = exp(— iif2^) and denote with Vie 
the coupling matrix acting on the Hilbert space Tii ® Tic- We therefore obtain: 

Vl,e{t) = Ui{ty ne(t)^ U2{ty [Vi^e ® I2 ] Ui{t) U^it) n2(t) 

= [ui{ty uo{ty i/i,c ui{t) uo{t) ] ® 12 • (40) 

This result can be interpreted as the "old" random coupling matrix in the interaction 
picture acting in the Hilbert space Hi 0He times the identity in 7^2- As a consequence: 

(^i,e(r) ViAr')) = Ci{t - r') 01^0 C,{r - t') (41) 

{Aj) = [ dr [ dr' [C,ir - r') Ci(r - r') ® I2 ^c + herm. conjg. ] .(42) 
Jo Jo 

In the case of (Aj) the calculation is a bit more involved: 

= Ce(r' -r)Y, Kj) e-(^-^')(^-^') ^^,,. ,„ (zn| . (43) 

ij,kn 

This result is quite similar to the expression in fl2Tl) for the same quantity in the non- 
structured case, where we had defined the operator-mapping Cc(t). One prominent 
feature of this mapping was the fact that it transformed any matrix into a diagonal one. 
In the present expression this occurs only partially, namely in Tii, while the density 
matrix remains in a sense untouched in 7^2- Thus we generalize Cc{t), in what follows 
Ci(t), to cases where its argument acts on a product space TCc = 'Hi0T-C2, in which case 
the operator mapping yields exactly what we obtained in Equation ( 143|) . 

Ci(i) : g, ^ CMit) = Y^ Kj) e-'(^-^^)(^-^') gl^^,^ {^n\ . (44) 

ij,kn 

This allows us to write: 

{Ai) = jj dr dr' Ce(r' - r) d [g,] (r - r') . (45) 

3.1. Master equation 

Following the same lines as the derivation of the master equation in the non-structured 
case, we find in the Fermi golden rule limit: 

{Aj) =miTHt g^ , {Ai) = r^ t li ® tri g^ . (46) 
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Therefore 

{Qc{t)) = &- A^r^t [mi ^c- li Otri^c] , (47) 

which may be continued in time, just as in the non-structured case, to give a closed 
differential equation for {gdt)): 

^ {0c{t)) = -mi th a' [ {g,{t)) - m^i Ui tri {Q,{t)) ] , (48) 

where (^c(O)) = ^c- 

Werner state solution Recently, [T7l [22]. it has been observed that the average purity 
and the average concurrence share a one-to-one correspondence during the decay of an 
initial Bell state. This correspondence is the same as for the one-parameter family of 
Werner states. Among other things, [T71 [22] consider the decoherence in a central two- 
qubit system coupled in the spectator configuration to a random matrix environment, 
and it is under these circumstances that the curious relationship can be observed. The 
master equation derived above resolves this puzzle, since, as we will show below, Werner 
states are solutions to that equation. The point is that the average purity (concurrence) 
agrees with the purity (concurrence) of the average density matrix in the limit of a 
large environment (A^ -^ oo). In the case of purity, this has been shown analytically 
within the linear response approximation and numerically in section HI In the case of 
concurrence this has been shown only numerically (in the same section H]). 

Let us consider the two-qubit system of Sec. HI We define a Werner state as a 
two-qubit mixed state which may be decomposed into a maximally entangled pure state 
and the identity matrix: 

g^(t) = ait)^ + Pit)\i;){i;\, (49) 

where ip is the maximally entangled initial state, and a{t) + (3{t) = 1. Maximal 
entanglement is equivalent to the condition 

tri|^)(^| = ^ ^ tr2|^)(V^| = ^. (50) 

We will show that gw{t) solves the master (H8l) for appropriately chosen coefficients 
a{t) and P{t). Note that (jUj) determines the time evolution of the two-qubits in the 
interaction picture. In order to obtain {gdt)) which describes the time evolution of the 
central system in the Schrodinger picture, we have to apply (TT2|) . If we have some single 
qubit dynamics, that may result in {gdt)) being no longer a Werner state. To complete 
the proof, we simply make the Ansatz: 

gdt) = Qw{t), gdO) = \^){ij\. (51) 

Substitution into the master equation fH5]) yields: 

« ^ + /? 1^) {^\ = -2th X' 13 (\^) {ij\ - y ) . (52) 
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Taking the trace on either side, we find: d + /? = 0, in agreement with the normahsation 
condition above. However, projecting either sides on |^)(^|, we obtain: 

| + /5 = -2r^A2^ => ^ = -2thX'P => /3(t)=e-2^«"'*. (53) 

As announced, this determines ait) and /3(t) in such a way that the Werner state Qwit) 
solves the master equation. 

3.2. Purity 

The hnear response expression ( ITSj) for purity equally holds for the present spectator 
model. Thus 



P(t) = tr (^,(t))2 ^ tigl - 2X^ ti[tT{{Aj)g,) - ti{{Aj)g,) ] , (54) 



where 



tT{{Aj)g,) = 11 drdr' C,{r - r') tr[Ci(r - r') ® 1^ gl] (55) 

ii{{Ai)g,) = jj drdr' C,{r' - r) tr[Ci[^,](r - r') g,] (56) 

If ^c represents a pure state, the purity P{t) in linear response approximation depends 
on the initial state only via tr2 ^c- This can be seen from: 

ti{{Aj)g,) = jj drdr' C,{t - r') tr[Ci(r - r') tra gl] , (57) 

which is true for any density matrix g^., and 

ii{{Ai)g,) = jj drdr' C,{r' - r) tr[Ci[tr2 g,]{r - r') ] , (58) 

which is only true if gc represents a pure state. 

4. Two-qubit central systems 

We shall proceed to analyse properties of the average density matrix for the particular 
case of a central system consisting of two qubits. Most of the analysis will be numerical, 
and thus we need a small central system for numerical expedience. A single qubit is of 
interest, but choosing this case would eliminate the option of considering entanglement 
within the central system. Anyway, the one qubit case will emerge as a special case, if 
the spectator is not entangled. 

For the two qubit case, we consider not only the purity of the density matrix p but 
also its von Neumann entropy S{p) = — trplogp. The latter quantity is preferred by 
many authors as a measure of decoherence. We shall quantify internal entanglement 
with concurrence, and study its properties in the spirit mentioned above. It is defined 
as C(p) = max{0, Ai — A2 — A3 — A4} where Ai are the eigenvalues of the matrix 
\/p{cry ® ay)p*{ay ® (Jy) lu non-increasing order, (*) denotes complex conjugation in 
the computational basis, and ay is a Pauli matrix 
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Figure 1. Difference between the entropy (squares) /concurrence (triangles) of the 
average density matrix and their value for each density matrix, averaged over the 
ensemble. The colours indicate different times and the lines correspond to a fit with 
fixed slope of — f. This shows that the average density matrix captures the physics 
of the ensemble for the quantities considered, including purity in the large A'' limit. 
Here, A — 0.03, and the average values of purity are P(rH/5) « 0.99, P(rH/2) « 0.95, 
Pirn) « 0.88 and P{2m) = 0.66. 



In figure [T] we show the differences between the ensemble averages of C and 5* as 
well as their value for the average density matrix p, for different times and stepwise 
increasing the size of the environment. As the reader can see this difference approaches 
zero with increasing environment size as 1/A^env Similar results were obtained for purity, 
but are not shown in a figure. They are consistent with the results found in section [2j 
The upshot is then, that we see all these very relevant quantities do not depend on the 
form in which the average is taken. Whether this consequence of a concentration of 
measure to some generalized delta function may be conjectured, but the demonstration 
thereof will be then interesting subject of future research. 

Concurrence and purity, for two non-interacting qubits coupled moderately to an 
environment, display a relation which corresponds well to that of a Werner state, while 
the members of the ensemble are not Werner states. It is thus reasonable to ask, whether 
the average density matrix corresponds to a Werner state. A 2 qubit Werner state, can 
be alternatively be defined as a density matrix with three degenerate eigenvalues, whose 
non-degenerate eigenstate is a maximally entangled state (i.e. a Bell state modulo local 
unitary transformations). We first examine under which circumstances we have a triple 
degeneracy. 

Numerically, we construct a finite ensemble of A^tot density matrices. We partition 
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Figure 2. Standard deviation of the tliree closest eigenvalues of the average density 
matrix as a function of the size of the ensemble (see text for details) . Different level 
splittings (differentiated by colours) and different times (differentiated by symbol type) 
are studied. When the coupled qubit has no internal Hamiltonian (blue symbols), 
the density matrix p has a triply degenerate eigenvalue in the large dimension limit. 
However when A > 0, this triple degeneracy is lifted, as can be seen from an asymptotic 
non-zero value. 



it in sets of equal size A^par, and evaluate the average density matrix p%\^^ for each of 

the Ntot/Nps,r members of the partition {i = 1, ■ ■ ■ , A^tot/A^par)- The properties of p^v ar 
are studied. In particular we consider the standard deviation awemer of its three closest 
eigenvalues. The average awemer for each of the Ntot/Npar sets available. This quantity 
is plotted in figure [2] for different level splitting A, [He = Aa^/^ in ([3])]. 

Indeed, if the coupled qubit has no internal dynamics, the average state, in the 
large dimension hmit, is triply degenerate. In fact awemer ^ as l/^/A^par- If the 
qubit has an internal Hamiltonian, the degeneracy is lifted: a systematic splitting of the 
eigenvalues (and hence a deviation from Werner states) happens on a time scale set by 
the mean level spacing of the environment. 

We also studied the concurrence of the eigenvector corresponding to the non- 
degenerate eigenvalue. It remains very close to one, independently of the concurrences 
of the other eigenvectors, which fluctuate without a clear pattern. We find the rather 
surprising result, that this property of the eigenfunction of the largest eigenvalue persists 
even in cases where the previous test shows, that we no longer deal with Werner states. 
Whether this property is characteristic of the average density matrix or whether it is 
common to many models will have to be analysed in future work. 
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5. Ensembles of density matrices and average density matrix 

The very spirit of the treatment of decoherence and entanglement by RMT leads to 
ensembles of density matrices: each member of the RMT ensemble (and, in our approach, 
also each initial condition in the environment) induces a density matrix pi via unitary 
dynamics followed by partial tracing. For that ensemble of density matrices we may 
consider the average density matrix p = p^. A difference of background will lead at this 
point to a natural inclination to two different approaches: Physicists with a background 
in decoherence, quantum optics and quantum information will most likely tend to use 
this average density matrix to calculate any quantity desired [say C(p)], while others 
with a background in random matrix theory will be inclined to calculate the ensemble 
averages of the quantities of interest [say C{pi)]. The difference between these two 
approaches disappears when one studies expected values of observables as tr Ap = tr Api 
or, as show in this paper, for some non-linear quantities in the limit of large dimension 
for the environment. Moreover, with quantitative entanglement witnesses [30], one 
can bound some of these non-linear quantities using (linear) observables. Yet the very 
concept of a relevant near environment, important for many of the new insights obtained 
form an RMT treatment, suggests to consider small environments as well. This leads 
us to discuss the appropriateness of either approach in typical applications. 

One must recall, that the determination of any of these quantities is not achievable 
as a single measurement. Rather we need an ensemble of "identical" systems. The 
hinge is what we mean by "identical" . If we have no control over the state or dynamics 
of the near environment, we must consider /(p) (/ being an expectation value, or any 
of the non-linear quantities discussed). Two notes are important at this point: First, 
the problem discussed here appears for any type of averaging, not only for averaging 
over time evolutions. Second, for the GUE the ensemble a state average is implicit in 
the ensemble average, while this is not the case for a GOE [31j, and thus has to be 
performed separately. However, if the near environment is known and its dynamics can 
be controlled, one should consider /(p). This can indeed be the case e.g. for an ion trap 
setup as the one discussed in the introduction, which might be used to experimentally 
verify our results. 

The differentiation between both approaches is also relevant within quantum 
information tasks. Consider for instance Shor's algorithm. The success rate has to 
be calculated from the average density matrix p resulting from many realizations of the 
same experiment. For the teleportation of an unknown state, however, it is important 
to know the success rate at each attempt, and thus the quality of the process should 
result from an average quality for individual density matrices corresponding to each try. 

At this point, to avoid confusion, we should recall that there is the option of 
obtaining a mixed density matrix without any entanglement by allowing probabilistic 
variations in the unitary time evolution of the central system itself in the absence of any 
significant coupling to the environment. Clearly we then have an ensemble of pure states, 
whose average naturally will not be pure. Yet errors due to variations of the unitary 
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evolution are usually considered as a loss of fidelity. The purity of the mixed density 
matrix would then measure an average quantity closely related to fidelity. This explains 
our emphasis throughout this paper on the fact, that we need a fixed evolution for the 
central system. It also shows, that the separation of external and internal perturbations 
for the study of the stability of quantum information tasks is somewhat artificial in a 
practical sense, though it is very useful for theoretical studies. 

The upshot of this discussion resides in the very fact, that for many quantum 
information tasks, we need repeated experiments. The decision, which average should 
be taken, is a subtle one. This emphasizes the importance of the result indicating, that 
the averages coincide in many cases. 

6. Conclusions 

We have presented a random matrix theory of decoherence and entanglement. The 
need for such a theory derives from two facts. On the one hand random matrix theory 
is known to provide a good generic description of properties of systems displaying, 
what is often known as quantum chaos or wave chaos. Indeed this statement is almost 
tautological as many authors nowadays omit any relation to classical chaos and use the 
relation to statistical properties of RMT as the definition of the latter. Fidelity decay 
is the other determining factor for the stability of quantum processes, in particular for 
those relating to quantum information. RMT has proven very successful in this field 
even in describing experiments [131 ISl [16] , and comparison with numerics for dynamical 
spin chains have also given encouraging results [321 ESI ES]. On the other hand the 
universal regime of exponential decay of coherence usually considered and derived in 
many ways (including RMT [Zl [8]) is basically founded in the typical situation of very 
long Heisenberg times in the environment. Intuitive arguments based on Fermi's Golden 
rule make this behaviour plausible almost independently of the model of the environment 
used, whenever decoherence occurs long before this time. Yet we argued, that situations 
with a near environment with fairly low level density and thus short Heisenberg times 
occur, and will become standard as quantum information systems with ever better 
isolation from the general environment are developed. It is under these circumstances, 
that RMT can provide the generic model to which the behaviour of specific systems 
should be compared. After describing the RMT model family, proposed to a large 
extent in earlier work, we proceed to analyse a point, which has set RMT models 
apart from other models of decoherence. The ensemble of evolution operators creates 
for the central system an ensemble of density matrices rather than a single density 
matrix. Consequently, properties such as entropy, purity and concurrence have so far 
been calculated as averages over that ensemble. However, we can also compute the 
average density matrix first, putting us on equal footing with other more conventional 
models. The various properties mentioned above are then determined from that single 
density matrix. This has also the great advantage to produce lower order quantities that 
have a fair chance to be calculated exactly using super- symmetric techniques. While we 
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have not yet achieved this goal, we have calculated the average density matrix in linear 
response approximation. This as well as numerics allowed us to compare purity, von 
Neumann entropy as well as concurrence of the average density matrix to the average of 
these quantities over the ensemble of density matrices. The central finding is, that for 
large environments at constant Heisenberg time (or mean level distance) the difference 
between the two approaches converges to zero as the inverse of the dimension of the 
relevant Hilbert space of the environment. This indicates, that for decoherence we 
can often use the average density matrix, and thus the RMT models are really on the 
same footing with usual descriptions. Yet we have to note, that, when describing the 
entanglement between two smaller systems, deviations are important and the question 
which average the behaviour of a single system should be compared to depends on the 
particular experimental situation. Consistently, if we specialize to a two-qubit central 
system, we find that average concurrence and the concurrence of the average density 
matrix will also approach the same limit for large environments. We found the more 
surprising fact, that at least for small decoherence the eigenfunction of the dominant 
eigenvalue of the average density matrix remains to very good approximation a Bell 
state, if the initial state was a Bell state. 

We have thus shown that purity and other quantities measuring entanglement yield 
the same result in the large environment limit, whether calculated from the average 
density matrix or as an average over the ensemble of density matrices. This was 
done analytically for purity in the linear response regime (section [2]) and numerically, 
for purity, von Neumann entropy and concurrence beyond the linear response regime 
(section Hj). 

Taking into consideration the convexity of this quantities, this might imply that 
in the large dimension limit of the environment the measure for the density matrices 
becomes similar to a Dirac delta in the sense that all or a large class of convex functions 
could be calculated directly from the average density matrix. However other results 
[T2] suggest that typical states are far from the average expected state. This apparent 
contradiction as well as the distribution of the density matrices as such shall be studied 
in a later paper. This work can also be readily extended considering more realistic 
RMT ensembles that RMT allows more realistic models then the classical ensembles. In 
particular the two-body random ensembles may play an important role particularly in 
their recent formulation for distinguishable spins |Mj. Also maps as generic models for 
gates are of possible interest, and we hope that we laid the foundation for more research 
in the domain of decoherence by "small" environments, where the Heisenberg time is 
shorter then or of the order of the decoherence time. 
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Appendix A. Average purity for the non-structured and the spectator 
model 

In this section, we review the results of [22] for the average purity instead of the purity 
of the average density matrix for the models treated in Sec. [2] and Sec. [31 We assume 
that the initial state is a product state of the form 

^0 = ^c® |V^e)(V^e| , (A.l) 

where the state of the central system ^c may be mixed, but the state of the environment 
\'ilje){ipe\ must be pure. The latter is a special requirement which allowed to map the 
spectator model onto the non-structured model in Ref. p2]. For the average purity it 
was found that 

{P{t)) = P{0)-\\Bj-B,) (A.2) 

5j = 4Rep[(J(t))f?o®f?o] (A.3) 

Bj = 2 {p[ (/(t) go I{t)) ® f?o ] - Rep[ (/(t) go ® I{t) go] + p[ {I{t) go ® go I{t) ]IA.4) 
p[A(g)B]=tTc{trcAtrcB) . (A.5) 

With Eq. dnD we obtain 

Bj = 4Re [ dr f dr' p[ (\/(r) V{t')) go® go] 
Jo Jo 

= 4Re [ dr [ dr' p[( Cc(r - r') ® Ce(r - r') ) f?o ® ^o] 
Jo Jo 

= 4 Re / dr / dr' tre( Cc{r - r') dr - r') g, g,) . (A.6) 

Jo Jo 

Since ge is a pure state: g"^ = g^, and because Cx(r — r') = C^ir' — r)*, we obtain: 

Bj = ARe f dr f dr' C,{r-T') Cc(r-r') =2 ff drdr' Cc(r-r') Cc(r-r') .(A.7) 

Jo Jo J Jo 

This is precisely the same quantity calculated in |22j whose final expression is given in 

Eqs. (A. 8) and (A. 16). It shows that the absolute values taken in [22], though correct, 

are not really necessary. 

To calculate the first term of Bj we note that 

(Vir) go V{r)) = \tj) e-'(^--^-)(--') {kl\ go \kl) {tj\ 

= Cc[f?c](r - r') Ce[f?e](r - r') . (A.8) 

Therefore, with Eq. ( I23l) . which holds equally well in the case of the environment: 

p[{I{t)goI{t))®go]= ff drdr'Cc(r'-r)tre[Ce[f?e](r-r') g,] 

drdr' C,{t' - r) 5e(r - r') , (A.9) 
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where S^ir — r') is defined in precise analogy to S^ij — r') in Eq. (1M|) . Similarly, we 
obtain for the following terms: 

p[ (/(t) go ® /(t) go) ] = II drdr' S,{t - r') S,{t - r') (A.IO) 



J JO 

V[ {lit) Qo ® Qo Ht)) ] = 11 drdr' S,{t - r') C,{r' - r) , (A.ll) 

where we had to assume again that gc represents a pure state. In summary, we obtain: 

{Pit)) = 1-2X^11 drdr' [Ce(r - r') C,{r - r') - S,{r - r') C,{r' - r) 

+ ^e(r - r') Se(r - r') - Ce(r' - r) S-^lr - r') ] , (A.12) 

where we have used the fact that the double integral over 5*0 (r — r') S'c(r — r') is 
automatically real. The whole result is obviously invariant with respect to an interchange 
of the subscripts e and c. This means that we could have defined p[A Cg> B] equally well 
by first tracing A and B over the environment and the resulting matrix product over 
the central system. Which is equivalent to the statement that the purity of the state 
of the environment (after tracing over the central system) is equal to the purity of the 
state of the central system (after tracing over the environment). 

If we compare the result for the average purity with Eq. (135|) for the purity of the 
average density matrix, we find that the difference is equal to 

(P(t))-P(t) = 2A2y^*drdr'5e(r-r')[S,(r-r')-Ce(r'-r)].(A.13) 

As discussed in [22], the function Se^r) shows a similar behaviour as Ce^r), except for 
an additional factor of order A^~^. This holds at least if the state of the environment is 
sufficiently delocalized in energy. For that reason, we expect the difference {P{t)) — P{t) 
to be of order A^~^. 
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